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ABSTRACT 

The single-degenerate model is currently a favourable progenitor model for Type la supernovae (SNe 
la). Recent investigations on the WD + He star channel of the single-degenerate model imply that 
this channel is noteworthy for producing SNe la. In this paper we studied SN la birthrates and delay 
times of this channel via a detailed binary population synthesis approach. We found that the Galactic 
SN la birthrate from the WD -I- He star channel is ~ 0.3 x 10~^ yr~^ according to our standard 
model, and that this channel can explain SNe la with short delay times 4.5 x 10^ — 1.4 x 10^ yr). 
Meanwhile, these WD -f He star systems may be related to the young supersoft X-ray sources prior 
to SN la explosions. 

Subject headings: binaries: close — stars: evolution — supernovae: general — white dwarfs 



1. INTRODUCTION 

Type la supernovae (SNe la) play an important role in 
astrophysics, especially in the study of cosmic evolution. 
They have been applied successfully in determining cos- 
mological parameters (e.g., Q, and A; Riess et al. 1998; 
Perlmutter et al. 1999). It is widely accepted that SNe 
la are thermonuclear explosions of carbon-oxygen (CO) 
white dwarfs (WDs) accreting matter from their com- 
panions (for a review see Nomoto et al. 1997). However, 
there is still no agreement on the nature of their progeni- 
tors (Hillebrandt & Niemeyer 2000; Ropke & Hillebrandt 
2005; Wang et al. 2008; Podsiadlowski et al. 2008), 
and this may raise doubts about the distance calibration 
which is purely empirical and based on the SN la sample 
of the low red-shift Universe. 

At present, various progenitor models of SNe la can 
be examined by comparing the distribution of the delay 
time (between the star formation and SN la explosion) 
expected from a progenitor channel with that of obser- 
vations (e.g., Chen & Li 2007; Xu & Li 2009; Lii et al. 
2009; Mannucci 2009; Schawinski 2009). Recently, there 
are three important observational results for SNe la, i.e., 
the strong enhancement of the SN la birthrate in radio- 
loud early- type galaxies, the strong dependence of the 
SN la birthrate on the colors of the host galaxies, and 
the evolution of the SN la birthrate with redshift (Delia 
Valle et al. 2005; Mannucci et al. 2005, 2006, 2008). The 
relation between SN la birthrate and radio power impli- 
cates the information on the time-scales of the order of 
10* yr, which corresponds to the radio activity lifetime; 
the strong dependence of the local birthrate on the col- 
ors of the host galaxies is related to the time-scales of 
the order of the galaxy color evolution (i.e., 0.5 — 1 Gyr); 
the evolution of the SN la birthrate with redshift is sen- 
sitive to the long time-scales (a few Gyr) (for details see 
Mannucci et al. 2006). According to the present observa- 
tional results, Mannucci et al. (2006) found that they can 
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be best matched by a bimodal delay time distribution, in 
which about half of the SNe la explode soon after star- 
burst, with a delay time less than ~ 10* yr, while those 
remaining have a much wider distribution, which can be 
well described by an exponential function with a decay 
time of about 3 Gyr (see also Mannucci 2008). Note that 
Scannapieco & Bildsten (2005) explored the two compo- 
nents of SN la birthrates and found that a young SN la 
population may be helpful to explain the Fe content of 
the intracluster medium in galaxy clusters. Moreover, 
by investigating the star formation history of 257 SN la 
host galaxies, Aubourg et al. (2008) recently found ev- 
idence of a short-lived population of SN la progenitors 
with lifetimes of less than 180 Myr. 

Over the last decades, two competing progenitor mod- 
els of SNe la were discussed frequently, i.e., the single- 
degenerate (SD) and double-degenerate (DD) models. Of 
these two progenitor models, the SD model (Whelan & 
Iben 1973; Nomoto et al. 1984; Fedorova et al. 2004; Han 
2008; Meng et al. 2009) is widely accepted at present. 
It is suggested that the DD model, which involves the 
merger of two CO WDs (Iben & Tutukov 1984; Webbink 
1984; Han 1998), likely leads to an accretion-induced col- 
lapse rather than a SN la (Nomoto & Iben 1985; Saio & 
Nomoto 1985; Timmes et al. 1994). For the SD model, 
the companion is probably a MS star or a slightly evolved 
subgiant star (WD -f MS channel), or a red-giant star 
(WD -f RG channel). However, these two SD channels 
did not predict such young SN la populations (Hachisu 
et al. 1996, 1999a, 1999b; Li & van den Heuvel 1997; 
Langer 2000; Han & Podsiadlowski 2004, 2006).'' 

Wang et al. (2009) recently studied a WD + He star 
channel for the SD model to produce SNe la. In the 
study they carried out detailed binary evolution calcula- 
tions of this channel for about 2600 close WD binaries 
with metallicity Z = 0.02, in which a CO WD accretes 
material from a He MS star or a He subgiant to increase 
its mass to the Chandrasekhar-mass limit. The study 

Note that Hachisu ot al. (2008) investigated new evolutionary 
models for SN la progenitors, introducing the mass-stripping effect 
on a MS or slightly evolved companion star by winds from a mass- 
accreting WD. The model can explain the presence of very young 
(< 10* yr) populations of SN la progenitors, but the model depends 
on the efficiency of the mass-stripping effect. 
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showed the SN la production regions in the (logP', Mj) 
plane (see Fig. 8 of Wang et al. 2009), where P' and M| 
are the orbital period and the mass of the He compan- 
ion star at the onset of the Roche lobe overflow (RLOF), 
respectively, and indicated that this channel is notewor- 
thy for producing SNc la. Because the WD + He star 
systems are from intermediate mass binary systems, this 
channel is likely to explain SNe la with short delay times. 
However, SN la birthrates and delay times through this 
channel are not well known from a viewpoint of the bi- 
nary population synthesis (BPS). 

The purpose of this paper is to study SN la birthrates 
for the WD + He star channel and to explore possible 
SN la progenitor systems with short delay times from this 
channel. In Section 2, we describe the BPS approach for 
the WD + He star channel. The simulation results of the 
BPS approach is shown in Section 3. Finally, discussion 
and conclusion are given in Section 4. 

2. BINARY POPULATION SYNTHESIS 

In order to investigate SN la birthrates and delay times 
for the WD + He star channel, wc have performed a se- 
ries of Monte Carlo simulations in the BPS study. In 
each simulation, by using the Hurley's rapid binary evo- 
lution code (Hurley et al. 2000, 2002), we have followed 
the evolution of 4 x 10^ sample binaries from the star for- 
mation to the formation of the WD + He star systems 
according to three evolutionary channels (Sect. 2.2). We 
assumed that, if the parameters of a CO WD + He star 
system at the onset of the RLOF are located in the SN 
la production regions in the (log P', M2) plane (Fig. 8 of 
Wang et al. 2009), a SN la is produced. Hereafter, we 
use the term primordial to represent the binaries before 
the formation of WD -|- He star systems. 

2.1. Common Envelope in Binary Evolution 

When the primordial primary (massive star) in a bi- 
nary system fills its Roche lobe, the primordial mass ratio 
(primary to secondary) is crucial for the mass transfer. If 
it is larger than a critical mass ratio, ^c, the mass trans- 
fer may be dynamically unstable and a common envelope 
(CE) forms (Paczyhski 1976). The mass ratio qc varies 
with the evolutionary state of the primordial primary at 
the onset of RLOF (Hjellming & Webbink 1987; Web- 
bink 1988; Han et al. 2002; Podsiadlowski et al. 2002). 
In this study we adopt qc — 4.0 when the primary is in 
the MS stage or Hertzsprung gap. This value is sup- 
ported by detailed binary evolution studies (Han et al. 
2000; Chen & Han 2002, 2003). If the primordial pri- 
mary is on the first giant branch (FGB) or asymptotic 
giant branch (AGB) stage, we use 

q, = [l.Q7-x + 2{^f]/2.13, (1) 

where Mf is the mass of the primordial primary, 

is the core mass of the primordial primary, and x = 
d In Ri I A In is the mass-radius exponent of the pri- 
mordial primary and varies with composition. If the mass 
donor stars (primaries) are naked He giants, q^ = 0.748 
based on equation (1) (see Hurley et al. 2002 for details). 

When a CE forms, the embedded in the CE is a 'new' 
binary consisting of the dense core of the primordial pri- 
mary and the primordial secondary. Owing to frictional 



drag within the envelope, the orbit of the 'new' binary 
decays and a large part of the orbital energy released in 
the spiral-in process is injected into the envelope (Livio 
& Soker 1988). The CE ejection is still an open problem. 
Here, we use the standard energy equations (Webbink 
1984) to calculate the output of the CE phase. The CE 
is ejected if 

V 2af 2ai / Ai?don 

(2) 

where A is a structure parameter that depends on the 
evolutionary stage of the donor, Afdon is the mass of the 
donor, Mace is the mass of the accretor, a is the orbital 
separation, Menv is the mass of the donor's envelope, 
Pdon is the radius of the donor, and the indices i and 
f denote the initial and final values, respectively. The 
right side of the equation represents the binding energy 
of the CE, the left side shows the difference between the 
final and initial orbital energy, and ace is the CE ejec- 
tion efficiency, i.e., the fraction of the released orbital 
energy used to eject the CE. For this prescription of the 
CE ejection, there are two highly uncertain parameters 
(i.e., A and Qco)- We usually set A to be 0.5 to constrain 
ace (de Kool 1990), although an exact calculation should 
take into account the issue that A depends on the stellar 
structure. In principle, wc expect < a^c < 1, but we of- 
ten find that Qco exceeds 1 for the purpose of explaining 
observed binaries. This may indicate that other energy 
sources may also contribute to the ejection of the enve- 
lope, e.g., the internal energy of the envelope (Han et al. 
1994, 1995; Podsiadlowski et al. 2003; Webbink 2008). 
As in previous studies, we combine ace and A into one 
free parameter aceA, and set it to be 0.5 and 1.5 (e.g., 
Lii et al. 2006). 

2.2. Evolutionary Channels to WD + He Star Systems 

According to the evolutionary phase of the primordial 
primary at the beginning of the first RLOF, there are 
three channels which can produce CO WD -h He star 
systems and then produce SNe la. 

(1) He star channel. The primordial primary first fills 
its Roche lobe when it is in the subgiant or RG stage 
(Case B mass transfer defined by Kippenhahn & Weigert 
1967). At the end of the RLOF, the primary becomes 
a He star and continues to evolve. After the exhaus- 
tion of central He, the He star which now contains a CO 
core may fill its Roche lobe again due to expansion of 
the He star itself, and transfer its remaining He-rich en- 
velope to the MS companion star, eventually leading to 
the formation of a CO WD -|- MS system. After that, 
the MS companion star continues to evolve and fills its 
Roche lobe in the subgiant or RG stage. A CE is possi- 
bly formed quickly because of dynamically unstable mass 
transfer. If the CE can be ejected, a close CO WD -|- 
He star system is then produced. The CO WD -|- He 
star system continues to evolve, and the He star may 
fill its Roche lobe again (due to orbit decay induced by 
the gravitational wave radiation or the expansion of the 
He star itself), and transfer some material onto the sur- 
face of the CO WD. The accreted He may be converted 
into C and O via He-shell burning, and the CO WD in- 
creases in mass and explodes as a SN la when its mass 
reaches the Chandrasekhar mass limit. For this channel, 
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SN la explosions occur for the ranges Mi^i 5.0—8.0 Mq, 



2,i 



2.0 



6.5 Mq and P' ~ 10 - 40 days, where Mi^i, 



M2,i and P' are the initial mass of the primary and the 
secondary at ZAMS, and the initial orbital period of a 
binary system. 

(2) EAGB channel. If the primordial primary is on the 
early AGB (EAGB, i.e.. He is exhausted in the center of 
the star while thermal pulses have not yet started), a 
CE will be formed because of dynamically unstable mass 
transfer. After the CE is ejected, the orbit decays and 
the primordial primary becomes a He RG. The He RG 
may fill its Roche lobe and start mass transfer, which 
is likely stable and leaves a CO WD + MS system. The 
following evolution of the CO WD + MS system is similar 
to that in the He star channel above, and may form a 
CO WD + He star system and finally produce a SN la. 
For this channel, SN la explosions occur for the ranges 
Mi,i ~ 6.0 - 6.5 Mq, M2,i ~ 5.5 - 6.0 Mq and P' ~ 
300- 1000 days. 

(3) TP AGB channel. The primordial primary fills its 
Roche lobe at the thermal pulsing AGB (TPAGB) stage, 
and the companion star evolves to a He-core burning 
stage. A CE is easily formed owing to dynamically un- 
stable mass transfer during the RLOF. After the CE ejec- 
tion, the primordial primary becomes a CO WD, then a 
CO WD + He star system is produced. The following 
evolution of the CO WD + He star system is similar to 
that in two channels above, i.e., a SN la may be pro- 
duced finally. For this channel, SN la explosions occur 
for the ranges Mia ~ 5.5 - 6.5 Mq, M2,i ~ 5.0 - 6.0 Mq 
and P' > 1000 days. 

2.3. Basic Parameters for Monte Carlo Simulations 

In the BPS study, the Monte Carlo simulation requires 
as input the initial mass function (IMF) of the primary, 
the mass-ratio distribution, the distribution of initial or- 
bital separations, the eccentricity distribution of binary 
orbit and the star formation rate (SFR). 

(1) The IMF of Miller & Scalo (1979, MS79) is adopted. 
The primordial primary is generated according to the 
formula of Eggleton et al. (1989) 



Mf = 



0.19X 



(1 - X)0-75 + 0.032(1 - X)0-25 ' 



(3) 



where X is a random number uniformly distributed in the 
range [0, 1] and Mf is the mass of the primordial primary, 
which ranges from 0.1 Mq to 100 Mq. The studies of the 
IMF by Kroupa et al. (1993) and ZoccaU et al. (2000) 
support this IMF. As an alternative IMF we also consider 
the IMF of Scalo (1986, S86) 



M^ = 0.3 



X 



1-X 



0.55 



(4) 



where the meanings of X and Mf are similar to that of 
equation (3). 

(2) The initial mass-ratio distribution of the binaries, 
q' , is quite uncertain for binary evolution. For simplicity, 
we take a constant mass-ratio distribution (Mazeh et al. 
1992; Goldberg & Mazeh 1994): 

n{q') = 1, < g' < 1, (5) 

where q' = /Mf . This constant mass-ratio distribu- 
tion is supported by the study of Shatsky & Tokovinin 



(2002). As alternatives we also consider a rising mass 
ratio distribution 



niq') = 2q', 0<q'<l, 



(6) 



and the case where both binary components are chosen 
randomly and independently from the same IMF (uncor- 
related). 

(3) We assume that all stars are members of binary 
systems and that the distribution of separations is con- 
stant in log a for wide binaries, where a. is separation and 
falls off smoothly at small separation: 



a ■ n{a) 



{a/aoY 



a < ao, 
gq < a < ai, 



(7) 



where asep « 0.07, aq = 10 Pq, ai = 5.75 x 10*^ Pq = 
0.13 pc and m « 1.2. This distribution implies that the 
numbers of wide binary systems per logarithmic interval 
are equal, and that about 50 percent of stellar systems 
have orbital periods less than 100 yr (Han et al. 1995). 

(4) A circular orbit is assumed for all binaries. The 
orbits of semidetached binaries are generally circularized 
by the tidal force on a timescale which is much smaller 
than the nuclear timescale. Moreover, a binary is ex- 
pected to become circularized during the RLOF. As an 
alternative, we also consider a uniform eccentricity dis- 
tribution in the range [0, 1]. 

(5) We simply assume a constant SFR over the last 
15Gyr or, alternatively, as a delta fimction, i.e., a single 
starburst. In the case of a constant SFR, we assume that 
a binary calibrated with its primary more massive than 
0.8 Mq is formed annually (see Iben & Tutukov 1984; 
Han et al. 1995; Hurley et al. 2002). Prom this calibra- 
tion, we can get SFR = 5MQyr~-'^ (see also Willcms & 
Kolb 2004). For the case of a single starburst, we assume 
a burst producing 10^^ Mq in stars. In fact, the SFR in 
a galaxy is neither a constant nor a delta function over 
the last 15 Gyr. A galaxy may have a complicated star 
formation history. We only choose these two extremes 
for a simplicity. A constant SFR is similar to the sit- 
uation of our Galaxy (Ymigclson & Livio 1998; Han & 
Podsiadlowski 2004), while a delta function to the situa- 
tion of elliptical galaxies or globular clusters. Under the 
assumption of the SFR as a delta function, one can ob- 
tain the delay time of SNe la from a progenitor channel, 
and then to compare with that of observations (e.g., Han 
& Podsiadlowski 2004). 

3. THE RESULTS OF BINARY POPULATION SYNTHESIS 
3.1. Birthrates of SNe la 

We performed six sets of simulations (see Table 1) with 

metallicity Z = 0.02 to systematically investigate Galac- 
tic birthrates of SNe la for the WD 4- He star channel, 
where set 1 is our standard model with the best choice of 
model parameters (e.g., Han et al. 2002, 2003, 2007). We 
vary the model parameters in the other sets to examine 
their influences on the final results. 

In Figure 1, we show Galactic birthrates of SNe la 
for the WD -|- He star channel by adopting Z = 0.02 
and SFR = 5MQyr~^. The simulation for our stan- 
dard model (set 1) gives Galactic SN la birthrate of 
~ 0.3 X 10~^ yr""'^, which is lower than that inferred 
observationally (i.e., 3 — 4 x 10^'^ yr"^; van den Bergh 
& Tammann 1991; Cappellaro & Turatto 1997). This 
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TABLE 1 

Galactic birthrates of SNe Ia for different simulation 

SETS, where set 1 IS OUR STANDARD MODEL. Qcc^ = CE 
EJEC;TI0N PARAMETER: n((j') = INITIAL MASS RATIO 
DISTRIBUTION; IMF = INITIAL MASS FUNCTION; eCC = 
ECCENTRICITY DISTRIBUTION OF BINARY ORBIT; V — GALACTIC 



BIRTHRATES OF SNe IA. 


Set 




n(q') 


IMF 


ccc 


ly {10--^ yr-^) 


1 


0.5 


Constant 


MS79 


Circular 


0.295 


2 


1.5 


Constant 


MS79 


Circular 


0.302 


3 


1.5 


Constant 


MS79 


Uniform 


0.284 


4 


1.5 


Constant 


S86 


Circular 


0.168 


5 


1.5 


Rising 


MS79 


Circular 


0.232 


6 


1.5 


Uncorrelated 


MS79 


Circular 


0.040 
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. . . set 3 
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Fig. 1. — The evolution of Galactic birthrates of SNe la for a 
constant star formation rate (Z = 0.02, SFR = 5M0yr~^). The 
key to the line-styles representing different sets is given in the upper 
left corner. The results of sets 2 and 3 almost coincide with that 
of set 1 . 

implies that the WD + He star channel is only a sub- 
class of SN la production, and there may be some other 
channels or mechanisms also contributing to SNe la, 
e.g., WD + MS channel, WD + RG channel or double- 
degenerate channel (see Meng et al. 2009 for details). 
Especially, as mentioned by Han & Podsiadlowski (2004), 
the WD + MS channel can give a Galactic birthrate of 
~ 0.6 — 0.8 X 10~^ yr~^, and is considered to be an im- 
portant channel to produce SNe la. 

According to the results of the six sets of Monte Carlo 
simulations, wc find that the BPS is sensitive to uncer- 
tainties in some input parameters, in particular the mass- 
ratio distribution. If we adopt a mass-ratio distribution 
for un-correlated component masses (set 6), the birthrate 
will decrease to be ^ 4 x 10^'^ yr^^, as most of the donors 
in the WD + He star channel are not very massive which 
has the consequence that WDs cannot accrete enough 
mass to reach the Chandrasekhar-mass limit. 

The SN la birthrate in galaxies is the convolution of 
the distribution of the delay times (DDT) with the star 
formation history (SFH) (e.g., Greggio et al. 2008): 

iy{t) = [ SFR{t - t')DDT{t')dt', (8) 

where the SFR is the star formation rate, and t' is the 
delay times of SNe la. Due to a constant SFR adopted 
in this paper, the SN la birthrate v{t) is only related to 




Iog(t/yr) 

Fig. 2. Similar to Fig. 1, but for a single star burst with a total 
mass of lO-'-'AfQ. The result of set 3 almost coincides with that of 
set 2. 

the DDT, which can be expressed by 

f 0, t<ti, 
DDT{t) = I DDT'{t), ti < i < ta, (9) 
i 0, t> t2, 

where ti and t2 arc the minimum and maximum delay 
times of SNe la, respectively, and the DDT' is the dis- 
tribution of the delay times between ti and ta- When t 
is larger than the t2, the equation (8) can be written as 

v{t) = SFR / DDT'{t')dt' = constant. (10) 

Therefore, the SN la birthrates shown in figure 1 seems 
to be so completely flat after the first rise. 

Figure 2 displays the evolution of SN la birthrates for 
a single starburst with a total mass of 10^^ Mq. In 
the figure we see that SN la explosions occur between 
~ 4.5 X 10'^ yr and ~ 1.4 x 10^ yr after the starburst, 
which can explain SNe la with short delay times (Scan- 
napieco & Bildsten 2005; Mannucci et al. 2006; Aubourg 
et al. 2008). The minimum delay time in the figure is 
mainly decided by the MS lifetime of a 8 Mq star (it is 
also the maximum mass for the progenitors of CO WDs). 
Moreover, after the primordial binary system evolves to 
a WD -I- He star system, the MS lifetime of the He com- 
panion star also contributes to the minimum time, but 
the time is short, e.g., the MS lifetime of a 1 Mq He star 
is only about 15Myr (Eggieton 2006). 

3.2. Distribution of Initial Parameters of WD + He 
Star Systems for SNe la 

Observationally, some WD + He star systems are pos- 
sible SN la progenitors (Wang et al. 2009). Further 
investigations are necessary for final confirmation of this 
(from both observations and theories). In this section, 
we will present some properties of initial WD -|- He star 
systems for SNe la according to our BPS approach, which 
may help to search for potential SN la progenitors. 

Figure 3 shows the distribution of the initial orbital 
periods of the WD -|- He star systems that ultimately 
produce SNe la with different ctceA. The simulation uses 
a metallicity Z = 0.02 and a constant initial mass-ratio 
distribution. The figure displays a result of the current 
epoch for a constant SFR. In the figure we can see that 
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Fig. 3. — The distribution of the initial orbital periods of the 
WD + He star systems which can ultimately produce SNc la for 
different OcoA. The simulation uses a metallicity Z = 0.02 and a 
constant initial mass-ratio distribution. The solid and the dotted 
histograms represent the cases with aceA = 0.5 (set 1) and OceA = 
f .5 (set 2), respectively. The number in every case is normalized 
to 1. 



1.0 



Fig. 4. — Similar to Fig. 3, but for the distribution of the initial 
masses of the CO WDs. 



there are obviously two peaks for each case. The left 
peak in these two cases results from the He star channel. 
Many of the SNe la in the right peak are also from the 
He star channel while others from the EAGB channel 
and the TPAGB channel (see Sect. 2.2). The He star 
channel has an important contribution to the formation 
of SNe la. This figure also shows that a high value of 
aceA leads to wider WD binaries, since a high value of 
accA is easier to eject the CE in the binary evolution. 

Figure 4 represents the distribution of the initial 
masses of the CO WDs. In the figure, a low value of 
aceA tends to have larger WD masses on average. The 
He star channel in Sect. 2.2, which allows stable RLOF 
to produce massive WDs (rather lead to dynamical mass 
transfer and a CE phase), is useful to understand this 
trend. According to our BPS simulations, we find that a 




1.6 2.0 
WHe(M®) 



Fig. 5. — Similar to Fig. 3, but for the distribution of the initial 

masses of the He donor stars. 

low value of aceA will produce more SNe la through the 
He star channel than other two channels, and then pro- 
duce more massive WDs on average. Figure 5 displays 
the distribution of the initial masses of the He donor 
stars. A low value of QcoA in the figure tends to have 
larger He star masses on average. This is also related to 
the stable RLOF, which leads to more massive compan- 
ion star resulting in the final larger He-core mass (He star 
mass). Moreover, a massive He donor star in the WD -|- 
He star channel will evolve more quickly and hence pro- 
duce a SN at an earlier time. 

When the WDs in Figure 4 increase their masses to the 
Chandrasekhar-mass limit, they will explode as SNe la. 
Meanwhile, the He donor stars in Figure 5, which afford 
material to the WDs though the RLOF, will loss a signif- 
icant amount of masses. We find that the He stars have 
masses ^ 0.6 — 1.7 Mq at the moment of SN explosions. 
Marietta et al. (2000) presented several high-resolution 
two-dimensional numerical simulations of the impacts of 
SN la explosions with companions. The impact makes 
the companion in the WD + MS channel lose a mass of 
0.15 - 0.17 Mq, but the impact in the WD + He star 
channel is still unknown. The companion in the WD -|- 
He star channel may lose more masses than that of the 
WD + MS channel. This is because the orbit separation 
at the moment of SN explosion from this channel is sig- 
nificantly less than that of the WD -|- MS channel, which 
may result in a much stronger impact to the companion. 
The surviving companion star from the WD -|- He star 
channel could be verified by future observations. 

4. DISCUSSION AND CONCLUSION 

Wang et al. (2009), based on equation (1) of Iben & 
Tutukov (1984), estimated the potential SN la birthrate 
through the WD -I- He star channel to be ^ 1.2 x 
10^"* yr^^ in the Galaxy. The birthrate from Wang et 
al. (2009) is higher than that in this paper, this is due to 
the fact that the long orbit period (i.e., > 1 day in Fig. 
8 of Wang et al. 2009) systems, which are considered 
to produce SNe la in eqiiation (1) of Iben & Tutukov 
(1984), do not contribute to SNe la in the simulation re- 
sults of this paper. In addition, Umeda et al. (1999a) 
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concluded that the upper hmit mass of CO cores born in 
binaries is about 1.07 Mq. If this value is adopted as the 
upper limit of the CO WD, the birthrate of SNe la from 
this channel will decrease to be ^ 0.2 x 10~^yr~^ in the 
Galaxy according to our standard model. 

In this paper we assume that all stars are in binaries 
and about 50 percent of stellar systems have orbital pe- 
riods less than 100 yr. In fact, it is known not to be the 
case, and the binary fractions may depend on metallicity, 
environment, spectral type, etc. If we adopt 40 percent of 
stellar systems have orbital periods below 100 yr by ad- 
justing the parameters in equation (7), we estimate that 
the birthrate of SNe la from this channel will decrease to 
be ~ 0.24 X lO^'^yr^^ for our standard model. 

SNe la from the WD + He star channel usually have 
massive CO WDs as their progenitors. Some previous 
studies showed that a massive CO WD leads to a lower 
C/0 ratio in the Chandrasckhar-mass WD, and thus a 
lower amount of ^^Ni synthesized in the thermonuclear 
explosion, which results in a lower luminosity of SNe 
la (Umeda et al. 1999b; Nomoto et al. 1999, 2003). 
However, brighter SNe la more frequently occur in ac- 
tive star formation galaxies (Hamuy et al. 1995, 1996), 
in which the young stellar population implies that these 
SNe la have short delay times (see also Aubourg et al. 
2008), i.e., the CO WD + He star channel might pro- 
duce brighter SNe la. Therefore, it is difficult to explain 
SN la diversity by using the C/0 ratio. Note that 3D 
simulations about SN la explosions by Ropke & Hille- 
brandt (2004) also indicated that different C/0 ratios 
have a negligible effect on the amount of ^^Ni produced. 
To understand the diversity of SN la explosions, the for- 
mation of brighter SNe la should be explored in future 
investigations. 

It is suggested that the WD -|- He star systems may 
appear as supersoft X-ray sources (SSSs) prior to SN la 
explosions (Iben & Tutukov 1994; Yoon & Langcr 2003: 
Wang et al. 2009). Recently, Di Stefano & Kong (2003) 
used a set of conservative criteria, applicable to Chandra 
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data, to identify luminous SSSs in four external galaxies 
(an elliptical galaxy, NGC 4967; two face-on spiral galax- 
ies, MlOl and M83; and an interacting galaxy, M51). 
They found that in every galaxy there are at least several 
hundred luminous SSSs with a luminosity of 10^ erg s~^ , 
and that in spiral galaxies MlOl, M83 and M51, the SSSs 
appear to be associated with spiral arms. This may indi- 
cate that some SSSs are young systems, possibly younger 
than 10* yr. Note that a WD -|- He star system has X-ray 
luminosity around 10^'' — 10^* erg s^^ when He burning 
is stable on the surface of the WD (Wang et al. 2009). 
Meanwhile, the distribution of SNe la with short delay 
times associated with galactic spiral arms (Barttmov et 
al. 1994; Delia Valle & Livio 1994). Therefore, we em- 
phasize that these WD -|- He star systems may be related 
to the young SSSs prior to SN la explosions. 

The most important conclusion of this study is that the 
WD -t- He star channel can explain SNc la with short de- 
lay times (~ 4.5 x 10^ — 1.4 x lO'^ yr), which is consistent 
with recent observational implications of young popula- 
tions of some SN la progenitors (Scannapieco & Bild- 
sten 2005; Mannucci et al. 2006; Aubourg et al. 2008). 
The young population of SNe la may have an effect on 
models of galactic chemical evolution, since they would 
return large amounts of iron to the interstellar medium 
much earlier than previously thought. It may also have 
an impact on cosmology, as they are used as cosmological 
distance indicators. 
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